Weighted scale-free network with self-organizing link weight dynamics 
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All crucial features of the recently observed real-world weighted networks are obtained in a model 



where the weight of a link is defined with a single non-linear parameter a as Wi 



(SiSj 



and Sj are the strengths of two end nodes of the link and a is a continuously tunable positive 
parameter. In addition the definition of strength as Si = SjWij results a self-organizing link weight 
dynamics leading to a self-consistent distribution of strengths and weights on the network. Using 
the Barabasi- Albert growth dynamics all exponents of the weighted networks which are continuously 
tunable with a are obtained. It is conjectured that the weight distribution should be similar in any 
scale-free network. 



PACS numbers: 89.20.Hh, 89.75.Hc, 89.75.Fb, 05.60.-k 
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Representing widely different large interacting systems 
as networks and identifying common generic characteris- 
tics in them made a significant advance in the study of 
complex systems in recent years |l|, |2|, y, til ■ However in 
most cases, link weights, representing the strength of tics 
between nodes, are not similar and may even be highly 
heterogeneous. E.g., number of passengers between two 
airports 0, 0] , strength of pair-interaction between two 
species in ecological system, number of coauthored pa- 
per of two authors, data traffic in a link of the Internet 
or volume of trade between two countries differ widely . 
Present availability of empirical data allows to study the 
weighted networks having weights associated with links. 
Weighted networks are usually described by the adja- 
cency matrix W, ij-th element Wij of which represents 
the link weight between i and j-th node. Recent anal- 
ysis of the World Airport Network (WAN) data shows 
that the link weights defined by the passenger traffic has 
power law distribution: P(w) ~ u> -7 ™ 0. 

Evolution of these weighted networks and the dynam- 
ics associated with it is a matter of general interest. In 
a traffic network like WAN when a new airport emerges 
with few links the passenger traffic associated with these 
new links affect the traffic of the other links also. In turn 
the effect of the other links also change the traffic of the 
new links and this mutual interaction takes its time to 
reach a steady state. Similar effects are there in many 
weighted networks. In this paper we try to explore this 
self-organisation of link weights and node strengths with 
the evolution of the network. 

A large number of models have been proposed to gener- 
ate different type of networks. Most well known among 
them is the Barabasi and Albert (BA) network which 
grows by sequential addition of new nodes using a lin- 
ear preferential attachment probability. These are un- 
weighted networks where each link has the same weight. 

Strength of a node is measured by the total amount 
of weight supported by the node: Sj = HjWij. The 
WAN data shows non-linear scaling between average 
strength of a node with its degree: s(k) ~ k 13 with 



[3 ~ 1.5. On the other hand average link weight again 
scales with the product of the degrees of the two end- 



nodes 



(kikjf, where 9 = 0.5 for WAN and 0. 



for E. Coli metabolic networks where link weights repre- 
sents the optimal metabolic fluxes pllii. Hof. 

Barrat, Barthelemy and Vespignani (BBV) proposed a 
model[lll| for the evolution of the weighted network simi- 
lar to the BA model. A new node n selects a node i of the 
network with a probability proportional to its strength: 
Hn—ti ~ Si. In addition, S amount of weight is distributed 
proportionally among the already existing links of the 
target node, hence with the topological evolution of the 
network link weights also evolve. BBV model produces 
power-law distribution of link weights, node degree and 
strength scales linearly and both exhibit power-law with 
same exponent. 

Even on an unweighted network one can assign link 
weights as ~ (kikj) e . This is a static assignment 
for a given structure. The idea is, large degree nodes 
process large amount of traffic and in the level of links, 
traffic through a link grows with the degree of its two 
end nodes. Using this idea we propose a model where 
the weight of a link depends jointly on the strengths of 
the end nodes which modifies the strengths of the nodes 
and which recursively modifies the link weights and so on 
until all link weights and strengths converges to a spe- 
cific set of fixed values. Apart from power law degree 
distribution we find power-law distributed probability of 
nodal strengths and link weights and the exponents de- 
cays with the model parameter a. The other crucial em- 
pirical observation for weighted networks (a) degree and 
strength of a node scales in a non-linear way and (b) a 
sub-linear variation of link weight with the product 
of the end-node degrees is also found to be present and 
the respective exponents and 9 varies with a. 

In our model, the growth of the network structure hav- 
ing weights associated with every link is executed by two 
mutually independent processes as follows: (i) The net- 
work is a BA network and is grown using the simple algo- 
rithm described in Starting from a small fully con- 
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FIG. 1: Scaling plot of probability distribution for (a) node 
strengths for TV = 2 12 ,2 14 ,2 16 for a = 0.40, (b) link weights 
for a = 0.40, for TV = 2 16 , 2 18 , 2 20 . 



nected m-clique, at every time step one link is selected 
with uniform probability and the new node is connected 
to one of the two end nodes with probability 1/2. (ii) 
After every node is added to the network, the weights of 
all links are re-determined by a self-consistent iterative 
process such that the weight w%j of the link between the 
i-th and the j-th nodes depends on the nodal strengths 



si and Si as: 



Si — Zjj Wij 



and 



(1) 



where a is a continuously tunable parameter. Indeed 
the dependence of link weight on the strengths of the 
end nodes has been observed in the International Trade 
Network with a ~ 1.33 |l3j . 

At every time step a new node i and a new link ij 
with unit weight is added to the network increasing the 
degree as well as the strength of the node j attached. 
As a result the weights of the other links meeting at j 
increases as well, which in turn enhance the strengths 
of the neighbors of j. In this way the effect of insert- 
ing a new node spreads throughout the network. One 
can imagine two widely apart time scales are involved in 
this process: the slow time scale characterize the network 
growth process whereas the fast time scale determines the 
re-organization of the link weights and strengths. 

Given a connected network, the distribution of weights 
over all links and the distribution of strengths over all 
nodes consistent with Eqn. (1) can be calculated in a self- 
consistent iterative process. Let s™ and wfj denote the n- 
th update of the strength of the i-th node and the weight 
of the ij-th link. In the self-consistent procedure, initially 
all links are assigned unit weights, ttfy = 1. The nodal 
strengths are then determined as: s° = YljW^. In the 
first stage of iteration process, the weights are determined 
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FIG. 2: (a) Variation of strength distribution exponent 7 S 
with parameter a and (b) variation of link weight distribu- 
tion exponent ^ w with a for TV = 2 16 in both cases. Circles 
represents the data points and the continuous lines are the 
theoretical prediction of eqn. (6) and the exponential fitting 
function, drawn for reference. 



as: 



and consequently the strengths as: 



Tjjwjj, and so on. More specifically 
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Under this successive iteration process, the distribu- 
tion of link weights as well as nodal strengths gradually 
converge to fixed values when further iterations do not 
change their values any more. In practice we associate a 
tolerance 6 = 1CP 6 for the maximal difference of weights 
in a single link between two successive iterations. As 
described above that in sf sums run up to the second 
neighbours which implies that estimation of s™ should 
have sums running up to the n-th neighbours to consider 
the contribution up to n-th shell neighbours. Since BA 
network has the small- world properties of small diame- 
ters (diameter T>(N) of the network grows logarithmically 
with system size TV) therefore from an arbitrary node the 
whole network is reached within a few steps. 

Let us first calculate how the mean nodal strength 
(s(k)} of the fc-degree nodes depends on the degree k 
as a function of the tuning parameter a. Let an arbi- 
trary node i have degree ki and the mean degree of its 
nearest neighbours be k nn . Following the self-consistent 
iterative procedure, the strength Sj of the node i after 
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FIG. 3: (a) Variation of strength degree correlation exponent 
(3 with a, open circles are data points and the solid line is 
the theoretical prediction of eqn.(5). (b) Plot of link weight 
degree correlation exponent 8 with a, N = 2 16 , solid line is 
the plot of the fitting function 6 — a + ba c . 



n-th iteration is: 

Si (n) = k 1+a{1+2a{1+2a{ n-th(l+a)))) 

ia(l+2a(l+2a(l+2a(...n-th)))) u 2(n-2)a" 
"•n \n 

Si {n) = k M ^k^k^~ 2)an (2) 

where n is an arbitrarily large integer. Now for a steady 
state value of Sj(n) both the series M(n) and D{n) must 
converge and reach a steady state value for any arbitrarily 
large n. The convergence condition gives: 

M(n) = l + a(l + 2a(l + 2a( n - th{\ + a)))) 

1 + (3) 



1 - 2a 



and 



D(n) = a(l + 2a(l + 2a(l + 2a(...n-tft)))) 
a 



1 - 2a 



(4) 



For finite steady state value of M(n) and D{n) both of 
the above equation leads to the condition: a < 0.50, be- 
yond which link weights and node strengths will diverge. 

Now for strength-degree dependence s oc k@ we assume 
k n has no significant dependence on k, which is true for 
BA topology, hence 



(3 « M(n) = 1 + 



1 - 2a 



(5) 



To find the exponent j s of the probability distribution of 
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FIG. 4: Plot of (a) Weighted average neighbour degree (k™ n ) 
with k for a =0.00, 0.10, 0.20, 0.25, 0.30, 0.35, 0.40, 0.45 
from bottom to top, serially (b) average disparity measure 
< k * Y(k) > with k for the same a values from bottom, to 
top, N = 2 14 . 



= If _\ ■ Which gives, 



7s = 1 + 



( 7fc - 1)(1 - 2a) 



1 



a 



(6) 



Probability distribution of strengths shows power law 
decay: P(s) ~ s _7s over several decades in Fig. 1(a). 
We find that a 'knee'-like step is developed in the dis- 
tribution near a ~ 0.25 and remains up to a ~ 0.45. 
With the increase of a from 0.25 the position of the knee 
shifted to wards the higher value of strength and become 
more prominent. In Fig. 1(a) we have shown the scaling 
plot of P(s) ~ s at a = 0.40, where the 'knee' is quite 
prominent. The scaling function : 



P(s) - N x g(s/N") 



(7) 



where A w 1.55 and fi w 0.93 are estimated giving "/ s = 
1.67. This value is consistent with the value found 
from slope measurement of the logP(s) vs log(s) plot 
(« 1.68). 

The probability distribution of the link weights again 



shows power-law variation: P(w) 



~ w 



-7™ 



For very 



node strengths P(s) 



~ s 



we use the general relation 



small values of a w 0.10 link weights are very small, 
even for a system of size N = 10 6 maximum value of u> 
is around 10 and power-law region is not well-defined. In 
Fig. 1(b) scaling plot of P(w) is shown for a = 0.40 and 
with N = 2 16 ,2 18 ,2 20 . Weight exponent j w , found from 
scaling analysis (= 1.45/0.93 w 1.56) is close to the value 
found from slope measurement of log P(w) vs log it; plot, 
(« 1.55). 

In fig. 2(a) variation of strength exponent 7 S with a 
is plotted. Open circles are experimental points and the 
solid line is the theoretical prediction. At larger values 
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of a maximum value of strength will increase and the 
probability of occurrence of large strength nodes should 
be increased. Drop in the exponent 7 S of the power- 
law decay represents this trend. The data points and the 
expression of 7 S both indicates the value of 7 S at a ~ 0.50 
as 1.0. From the measurement of the exponent (3 also it 
was found that the value of j3 is asymptotically large as 
a — ► 0.50. In fig. 2(b) variation of link weight distribution 
exponent j w with a is shown. The initial high j w value 
shows that indeed for small a the distribution decays very 
fast. Data points are compared with a fitting function: 
7m = a. exp(— b.a), a = 7.98,6 = 4.16. This exponential 
decay of j w again signifies that link weights grow fast 
with a in general but specially the growth of large links 
are much faster. 

Fig. 3(a) shows strength-degree exponent (3 increases 
with a. While (3=1 implies that link weights are uncor- 
rected, for (3 > 1 large weights prefer to join the large 
degree nodes. Increase in the value of (3 means more 
aggregation of large weight links on large degree nodes. 
That signifies that a systematically increases the prefer- 
ence of large weights to converge on high degree hubs. 
The other factor is that with the increase of a degree 
values of the nodes remains unchanged but links of very 
large weights and hence very large node strengths appear 
and contributes to f3. The plot show a steep rise in in 
the region a > 0.40. In fig. 3(b) link weight-end de- 
gree exponent 9 is plotted with a. Exponent 9 found to 
grow with a. It is expected here as with a link weights 
increases but the degree distribution and hence the de- 
gree of any two end nodes remains unchanged. Plot is 
fitted with a function of the form: 9 = a + ba c , with 
a = 0.06,6 = 5.23, c = 2.07 where a « 0.00 is more 
appropriate as it corresponds to a = 0.00 value. 

The degree-degree correlation (k nn ) is 0: 

k n n,i — "T'^-'j^j \P) 

where the sum j runs over the nearest neighbours of i. 
Now when averaged over the nodes according to degree k 
we find k nn {k) which is related to the conditional proba- 
bility P(kx\k) as that a node of degree k is connected to 
a node of degree k±: 

knn(k) = -^^k,=kknn,i = Sfej klP(k% \k) (9) 

If there is no degree-degree correlation P(ki\k) is a func- 
tion of fci only and k nn (k) is a constant. If k nn (k) 
increases (decreases) with k, the network is said to 
be 'assortative' ('disassortative'). In weighted networks 
weighted average nearest neighbour degree is defined as 
follows: 

Cm = ->v"'<^ ( 10 ) 

and similarly the k^ n (k). 



It is known that in BA model there is no significant 
degree correlation. This has been observed in our cal- 
culation also that in Fig. 4(a) for a = 0, the curve is 
nearly horizontal for a large portion of the whole range of 
variation of the degree k. Therefore the underlying BA 
network structure cannot be responsible for additional 
correlation in k™ n (k) in Fig. 4(a) for a > and it must 
be due to the self-consistent procedure with which the 
links are updated. 

Our understanding is that the assignment of link 
weights as w^j oc (siSj) a induces this assortative na- 
ture. As initially Sj = fcj, weights of the links between 
large degree nodes become large and that produces the 
assortative nature in the weighted degree correlation as 

h w — iV .a/), h - 
h nn,i ~ Si 3 Wl 3 3 • 

The nodal disparity or inhomogeneity is Yi = 
Yij{wij j Si) 2 . When all the links contributed equally 
Yi ~ 1/fcj and when one link gives the major contribu- 
tion and others are not significant then Yi ~ 1. Network 
disparity Y{k) is the nodal disparity averaged over all 
nodes of degree k, shown in Fig. 4(b). For small values 
of a the plot is almost horizontal and the value is 1, in- 
dicating homogeneity of links for different node degrees. 
But for a > 0.35 there is inhomogeneity within a certain 
region of k value. First it starts growing then decreases 
to 1 again for hub nodes. This signifies that for a broad 
region of k values nodes have highly heterogeneous link 
weights but the links of the hub nodes are in general have 
very large weights and not heterogeneous in nature, that 
again is the effect of large link weights aggregation on 
hubs. 

To summarize, we have shown that starting from an 
unweighted network one can have non-trivial distribution 
of link weights as well as nodal strengths using only the 
definition of strength as well as the relation describing the 
non-linear dependence of link weights on the strengths 
of end nodes. The distribution of weights defines the 
strengths which in turn modifies the link weights - as 
a result link weights and nodal strengths self-organise 
themselves to a steady-state distribution through a re- 
cursive process. Study of this process on the BA network 
shows rich properties like power-law distributed probabil- 
ity distribution of nodal strengths as well as link weights, 
with associated exponents varying continuously with the 
model parameter a. The model exhibits all non-trivial 
features of real- world weighted networks, e.g., non- linear 
strength-degree relation, correlation between link weights 
and the product of degrees of the end-nodes etc. We con- 
jecture that similar non-trivial weight and strength dis- 
tributions can be obtained in general for any scale-free 
network and therefore may be useful to study different 
real- world weighted networks. 
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